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Abstract 

Noncommutative U(M) gauge theories at different A/ - may be often thought of as 
different sectors of a single theory: the U(l) theory possesses a sequence of vacua la- 
beled by an integer parameter M, and the theory in the vicinity of the Af-ih vacuum 
coincides with the U{N) noncommutative gauge theory. We construct noncommuta- 
tive domain walls on fuzzy cylinder, separating vacua with different gauge theories. 
These domain walls are solutions of BPS equations in gauge theory with an extra term 
stabilizing the radius of the cylinder. We study properties of the domain walls using 
adjoint scalar and fundamental fermion fields as probes. We show that the regions on 
different sides of the wall are not disjoint even in the low energy regime — there are 
^r) ' modes penetrating from one region to the other. We find that the wall supports a chiral 

fermion zero mode. Also, we study non-BPS solution representing a wall and an anti- 
wall, and show that this solution is unstable. We suggest that the domain walls emerge 
as solutions of matrix model in large class of pp-wave backgrounds with inhomogeneous 
field strength. In the M-theory language, the domain walls have an interpretation of a 
stack of branes of fingerstall shape inserted into a stack of cylindrical branes. 

1 Introduction 



Oh 
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Recently, field theories on noncommutative (NC) spaces attracted considerable interest (see, 
e.g., Refs. [1, 2, 3, 4] for recent reviews and references to earlier works). One of the reasons 
for this interest is the fact that NC gauge theories are non-local and that the group of gauge 
transformations in NC theories contains some of the diffeomorphisms 1 . Both non-locality 
and invariance under diffeomorphisms are expected to hold in the theory of quantum gravity, 
so one may hope to gain some insight into techniques appropriate in quantum gravity, using 
NC theories as a toy model (recall that NC theories are simpler in the sense that non-locality 
is present there already at the classical level). This hope is further supported by the fact 
that NC gauge theories emerge as effective descriptions of string theory in a certain limit 
[5], and are inherent in the matrix approach to M-theory. 

1 Strictly speaking, the latter property holds only for fields in the adjoint representation of gauge group. 
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Among the consequences of non-locality and invariance under coordinate transformations 
is the background dependence of both the matter content and space-time interpretation of 
NC gauge theory. One illustration of this background dependence is that in many cases 
NC gauge theories with different gauge groups U (M) emerge as different sectors of a single 
theory [6, 7, 8, 9]. For instance, it was pointed out in Ref. [6] that U(l) gauge theory on the 
NC plane possesses a sequence of vacua labeled by a natural number Af with the following 
peculiar properties: 

i) Every vacuum with Af > 1 is a highly non-local field configuration from the point of view 
of the trivial (A~ = 1) vacuum; 

ii) Perturbation theory in the vicinity of the ./V-th vacuum is equivalent to perturbation 
theory of the U (M) NC gauge theory above its trivial vacuum; 

iii) The fact that there are different gauge theories in different vacua cannot be understood 
as Higgs mechanism. Namely, the action in the vicinity of the Af-th vacuum contains U (Af) 
gauge fields only, with no extra massive vector bosons or Higgs fields. 

The precise physical meaning of this phenomenon is not clear yet. In the case of NC 
plane, it was argued [6] that Af is a superselection parameter, implying that different sectors 
are completely disconnected from each other. 

One way to understand the physical meaning of the parameter A" is to study whether 
there exist domain walls separating vacua with different values of A". These walls, if any, are 
expected to exhibit interesting physical properties, since in the low-energy (commutative) 
limit they would serve as boundaries between regions with different gauge theories inside. 
On the other hand, the absence of the domain walls would imply that sectors with different 
Af are disconnected, and A" is indeed just a superselection parameter. 

In this paper we study the problem of existence of the domain walls in the context of 
gauge theory on a fuzzy cylinder introduced in Ref. [8] (see also Ref. [10]). We choose 
cylindrical geometry because in this case a domain wall normal to the axis of the cylinder 
has finite energy. Fuzzy cylinder can be thought of as a lattice version of the continuous 
NC cylinder studied in Ref. [9]. It is a well-defined NC space which reduces to ordinary 
cylinder in the commutative limit. An advantadge of the fuzzy cylinder setup compared to 
the continuous NC cylinder is that the intriguing phenomenon described above (existence of 
the sequence of vacua corresponding to different gauge theories) is more transparent in the 
former case. We elaborate on the relation between the fuzzy and NC cylinders in section 2. 

We find that though configurations of finite energy interpolating between different gauge 
vacua exist already in pure Yang-Mills theory on the fuzzy cylinder (these configurations 
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were first constructed in Ref. [11] and were interpreted there as D-brane junctions), they are 
unstable and tend to dissolve classically 2 . However, once a gauge-invariant term, stabilizing 
the radius of the cylinder, is added to the action, the domain walls become stable solutions. 
They saturate the BPS bound emerging due to the existence of a non-trivial topological 
charge. We construct explicitly a family of solutions of the BPS equations, that describe 
domain walls between U{M\) and U(Nq) vacua with arbitrary M\ and A/2. 

We suggest also a matrix model interpretation of these domain walls. Namely, we argue 
that a term needed to stabilize the domain wall is generated in a large class of curved 
gravitational pp-wave backgrounds with inhomogeneous three-form field strength. In the 
M-theory language, domain wall solutions have an interpretation of a stack of branes of 
fingerstall shape, inserted into a stack of cylindrical branes. 

An interesting property of the domain walls we construct is that the regions with different 
gauge theories they separate are not disjoint even in the low energy (long wavelength) regime. 
We illustrate this point by studying the properties of adjoint scalar and fundamental fermion 
fields in the background of the simplest domain wall between U(l) and U(2) gauge vacua. 
We show that there are long wavelength modes of these fields, which penetrate from the 
region with U(2) gauge theory to that with U(l) gauge theory. It is worth stressing that 
these modes do not belong to the diagonal part of the U(2) group. Other modes, on the 
contrary, experience total reflection from the wall. In addition, we find that the wall localizes 
a zero fermionic mode, whose profile can be used as a probe of the shape of the wall. 

Finally, we study a wall-antiwall system which separates the cylinder into three regions 
with U(l), U(2) and again U(l) gauge theories. An unusual property of this system (which 
is not BPS) is the absence of attraction between the wall and antiwall in the sense that there 
is a one-parameter family of solutions with equal energies and different distances between 
the wall and antiwall. However, a tachyonic mode is present in the spectrum of excitations of 
the wall-antiwall system. The absolute value of its mass squared is exponentially small when 
the distance between the walls is large. Thus, the wall-antiwall system is almost stable, 
when the walls are well separated. 

This paper is organized as follows. In section 2 we introduce the algebra of functions on 
the fuzzy cylinder and describe the gauge theory. In section 3 we discuss the topology of 
domain wall configuration, derive the BPS bound for its energy and solve the correspond- 
ing BPS equations. In section 4 we consider scalar and fermion fields in the domain wall 
background. Wall-antiwall system is studied in section 5. In section 6 we discuss the rela- 
2 In Rcf. [11] it was suggested that these configurations may be stabilized by quantum effects. 
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tion between our model and matrix theory in curved background. The concluding section 7 
contains brief summary and discussion of our results. 

2 Fuzzy cylinder 

2.1 Algebra of functions on the fuzzy cylinder 

To introduce algebra Ac of functions on the NC cylinder and algebra Af of functions on 
the fuzzy cylinder, it is convenient to start with the "master" algebra A generated by three 
elements x, y, z, obeying the following commutation relations 

[z, x] = ily , [z, y] = -iZx , [x, y] = 

where the parameter / is the scale of noncommutativity and operators x, y and z may be 
thought of as the coordinates of a three-dimensional NC space where NC and fuzzy cylinders 
are embedded into. It is convenient to introduce linear combinations 

x+ = x + iy , x = x - iy . 

Commutation relations for these elements have the form 

[z, x+] = Zx+ , [z, x_] = -Zx_ , [x+, x_] = . 

Irreducible representations of the algebra A are labeled by the eigenvalues of the two central 
elements 

Ti = x+x = x 2 + y 2 , T 2 = e 2wiz/l . 

The eigenvalue of the first central element Ti is natural to interprete as radius squared of 
a cylinder. Thus, to introduce algebra Ac of functions on the NC cylinder of radius p, one 

sets 

Ti = p 2 . 

More formally, algebra Ac can be defined as a factor-algebra 

Ac = A/{Ti - p 2 } , 

where {Ti — p 2 } is a subalgebra of algebra A generated by the element (Ti — p 2 ). For recent 
study of gauge theory on the NC cylinder see Ref. [9]. 
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On the other hand, by fixing the value of the second central element 

one obtains a collection of planes parallel to (x, y)-plane with z-coordinate equal to (z + nl), 
n G Z. Algebra Af of the fuzzy cylinder (see, e.g., Refs. [8, 10]) is obtained by fixing the 
values of both central elements 7\ and T 2 . In other words, one defines algebra Af as the 
following factor-algebra, 

A F = A C /{T 2 - 1} = A/{T X - p\ T 2 - 1} , 

where without loss of generality we set z = 0. This algebra may be realized as the algebra 
of operators acting in a Hilbert space H with basis vectors \n),n G Z, 

oo oo oo 

z = I ^ n|n)(n| , x + = p ^ |n + , x_ = p ^ \n — l)(n\ . (1) 

n=— oo n=— oc n=— oo 

It is clear from the above discussion that fuzzy cylinder can be thought of as a cylindrical 
semi-lattice with continuos coordinate 6 defined by 

x ± = pe ±w 

and discrete coordinate z with spacing /. 

To make this picture more transparent it is instructive to introduce symbols of the op- 
erators in the Hilbert space H. We will use the symmetric ordering which maps functional 
exponents into operator exponents 



e i(kz+N0) _^ e i(kz+N6) _ e -ikNl/2 e ikz 



X , 



Then it is straightforward to obtain the following relation between an arbitrary operator 
/ = Y.fnm\n)(m\ and its symbol f(z,9) 

oo r- ji 

/>,*)= * E fnrn e^~^)^-^ (2) 



n,m=— oo l 



Note that the integration over k has finite range k e (— tt/1, 7r/Z) . Consequently, the symbol 
f(z, 6) is uniquely determined by its values on the cylindrical semi-lattice with lattice points 

z = nl , n G Z . 
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Equivalently, one may view the map (2) as a correspondence between operators in H and 
functions on the cylinder, whose Fourier components along z-coordinate are cut off at the 
scale ir /I. Clearly, a symbol f(z,9) considered as a function of the continuous variable z 
contains the same amount of information as its values on the lattice f(ln,9), but many 
formulae simplify when written in terms of functions of the continuous variable. 

The map (2) implies the following relation between trace of the operator and integral of 
the symbol, 

2nplTrf = I pdOdzf = plJ2 [ d9f(ln, 9) . (3) 

Let us now define derivatives of functions on the fuzzy cylinder. It is straightforward to check 
that differentiation of symbols with respect to ^-coordinate translates into the following inner 
derivation in the operator language, 

tf«{E7i-g. (4) 

This relation has the same form as in the case of NC cylinder (see, e.g. Ref. [9]). However, it 
is impossible to define a derivative along ^-direction. Indeed, in the language of symbols we 
consider functions on the lattice in z direction, so it is natural to expect that some discretized 
version of the derivative in z direction emerges. In the operator language, a naive attempt 
to define z-derivative would contradict the Leibnitz rule because of the constraint T 2 = 1. 

On the other hand, from the algebraic point of view it is natural to consider on equal 
footing the following three inner derivations of the algebra Af 

fc/ = }[x,/], d 2 f = -j[y,f], d 3 f = j[z,f]. 

As pointed out above, the last derivation 83 corresponds to differentiation with respect to 9 
in terms of symbols. Furthemore, it is straightforward to check that 

dj= y[x,/] =yd z f (5a) 

= -jfof\ = xdj , (5b) 

where finite-difference derivative d z is defined as 

f(z + l/2)-f(z-l/2) 



dj 



I 



Note that this derivative has a simple form when written in terms of symbol f(z) considered 
as a function of continuous variable z. However, written as a lattice derivative this operator 
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has the following highly non-local form 



/4 

This property illustrates the fact that the formulation in terms of functions f(z) of continuous 
variable is often more convenient than the lattice formulation. 

2.2 Scalar field on the fuzzy cylinder 

To get accustomed to physics of fuzzy cylinder, let us consider a free scalar field theory on 
it. The action has the following form 



S = J dt 2ixpl Tri (V o 0) 2 - -^((5i0) 2 + (d 2 0) 2 + (9 3 0) 2 ) - m 2 2 ) 



(6) 



The spectrum of this theory can be determined in two different ways. First, one can rewrite 
the action (6) in terms of the symbol of the operator <fi. Using Eqs. (3) - (5) one has 

S = J dt pdQdz 1 - (Vo0) 2 - ^(<%0) 2 - (40) 2 - m 2 ~^j . (7) 



Solutions of the field equations following from the action (7) have the form of waves 

7T 7T 

7' 7 



oc e - Mkz+im , ke' " " 



with the dispersion relation 



2 N 2 (2 . klV 2 
- 2 = -2-+ ysm- + m 2 . (8) 



In the long wavelength limit Ic < 1/1 we recover the usual dispersion relation for scalar waves 
on the cylinder. 

It is instructive to obtain the dispersion relation (8) in the operator approach. Variation 
of the action (6) yields the following field equation 

= ^2 (j Z ' Ml] + ^ X +> [ X -'^] + ^ X "> [ X +><^]) +™V • 

In components nm of the operator = (f> nm \n)(m\, this equation takes the form 

1 / ] 

— (n-m) (p nm + y 



-dl<p nm = \ (« ~ ™) 2 0nm + 4(20nm ~ 0ra+l,m+l - 0n-l,m-l) + ™Vnm • (9) 
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The system of equations (9) decomposes into independent recursion equations along diago- 
nals of the matrix 0. This corresponds to the Kaluza-Klein decomposition over the compact 
variable 6. Let us consider equation (9) along the iVth diagonal, n — m = N, describing the 
Nth KK-mode, 

-<9q0„ = -^N 2 <p n + ^(20„ - (f) n+1 - n _x) + m 2 (j) n . 

where we have set 4> n = 4> ntn+ N- We immediately find that the Ansatz <p n = oe -^'+ lfcfa [ s 
consistent with this equation, yielding the dispersion relation (8). 

2.3 Gauge theory on the fuzzy cylinder 

Let us now describe U(Af) gauge theory on the fuzzy cylinder. To this end we consider a 
field ip transforming under the fundamental representation of U(Af). This field belongs to 
the direct sum of Af copies of the algebra Af 

iPeC M (g)A F = ^A F . 

Gauge transformations are defined as follows, 

where U is a unitary operator acting in C^" <8> Af (equipped with scalar product (^1,^2) = 
Tr(ip^ip2))- Covariant connection is introduced in the following way 

V(9x)^ = y(X^-^x) 

v(W = yW-v>z), 

where covariant coordinates X, Y, Z are Hermitian Af x Af matrices with entries in Af- 
They transform under the adjoint representation of the gauge group, X — > U~KU + , etc. The 
covariant strength tensor is defined in the usual way, 3 

Fa = [y (0Q, v(a,-)l - v([ft, dj\) . (10) 

3 The second term in the r.h.s. of Eq. (10) may appear somewhat unusual. In fact, it is present in the 
conventional field theory as well. In commutative field theories this term ensures that is a tensor function 
(not a differential operator) when vector fields 9, and dj do not commute. 
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Thus, we obtain 




(lid) 



(lib) 



(11a) 



(11c) 



Now, it is straightforward to write down the Yang-Mills action, 



(12) 



where summation is assumed over indices i,j = 1, 2, 3. To figure out the commutative limit 
of the theory, it is convenient to decompose the covariant coordinates in the following way, 



Substituting expressions (13) into Eqs. (11) and taking the limit I — > 0, one obtains 



F a/3 = d a A p - d p A a - i[A a ,A p ] , D a A p — d a A p — i[A a , A p ] , a,(3 = t,9,z. 

Thus, in the commutative limit one obtains the usual Yang-Mills theory on the cylinder, 
coupled to an adjoint scalar field A p . The presence of such a scalar is a consequence of 
the fact that the algebra of the fuzzy cylinder Af possesses three independent derivations 
instead of two, and hence there are three gauge fields in the NC gauge theory. 



X = x - -(yA z + A z y) + -(xA p + A p x) 
Y = y + 1 -{*A Z + A z x) + l -{yA p + A p y) 

Z = Z - lAg . 



(13b) 
(13c) 



(13a) 



F 12 = -ip 2 D z A p 



F\z = ip(sm0Fg z - cos 9D g A p ) 
F 23 = ip(cos 9F ez + sin 6D e A p ) 
F m = ip(- sin 9F tz + cos 9D t A p ) 
F 02 = ip(cos9F tz + sm9D t A p ) 



F03 = —iF t e , 



where 
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Now, let us note the following peculiar property of NC gauge theory on the fuzzy cylin- 
der 4 . Consider U(l) gauge theory. Then for any natural Af, there exists a vacuum in this 
theory, such that the theory above this vacuum is identical to U (Af) gauge theory above its 
trivial vacuum. Namely, the vacuum corresponding to U (Af) theory has the form 

oo M 

Z = l ^ n^2\nAf + a)(nAf + a\ (14a) 

n=— oo a=l 
oo M 

X + = pJ2 ^|(^+l)AT + a)(nAT + a| (14b) 

n=— oo a=l 

A = . (14c) 

It is straightforward to check, that the energy of configuration (14) is equal to zero, and thus 
it indeed describes a vacuum in the U(l) theory. Equivalence between this vacuum and the 
trivial vacuum in U (A/*) theory can be established using an isomorphism S : <E> H — > H , 
defined on the basis vectors as 

S 

| a) <S> \n) I— > \nAf + a) 

which maps operators (14a), (14b) into the following operators acting in a direct sum of M 
Hilbert spaces 

/ oo \ 

(15a) 

1 , (15b) 

where 1 stands for the unit Af x Af matrix. Field configuration given by Eqs. (15) describes 

a trivial vacuum in the U(N) theory. Using the same trick for the field configurations 

describing fluctuations above the vacuum given by Eqs. (14), one observes that the action 

governing these fluctuations is equivalent to the action describing fluctuations of the U (AT) 

theory in the vicinity of the trivial vacuum. 

Thus, the U(l) gauge theory on the fuzzy cylinder has an infinite set of vacua labeled by 

Af = 1, 2, . . . It is natural to wonder whether these vacua correspond to different superselec- 

tion sectors, or are different phases of one and the same theory. In particular, one may ask 

4 This property is generic for gauge theories on a non-compact NC manifold. In the case of NC plane it 
was discussed in Rcf. [6] , for NC cylinder in Ref. [9] ; for earlier discussions of this property for fuzzy cylinder 
see Refs. [8, 11]. 
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whether there exist solutions (domain walls) interpolating between different vacua. In the 
next section we answer affirmatively to this question and explicitly find such domain wall 
solutions. 

3 Domain wall 

In what follows it is convenient to use the dimensionless variables 

Z=\z, X ± = -X ± . (16) 

We search for a static solution of U(l) gauge theory on the fuzzy cylinder in the form of a 
domain wall, separating vacua corresponding to different gauge theories. We will work in the 
gauge A = 0. For concreteness we first concentrate on the case of a domain wall between 
U(l) gauge theory in the region z < and U(2) gauge theory in the region z > 0. Then the 
asymptotics of the domain wall are given by Eqs. (1) and (14), respectively, 

Z = n\n) (n\ , X + = \n + 1) (n\ , n — > — oo (17a) 
Z = ^n(|2n-l)(2n-l| + ]2n)(2n|) , X + = ^ \n + 2)(n\ , n -> +oo . (17b) 

Let us note that field configurations of finite energy, possessing asymptotics (17) do exist 
[11]. A simple example is 

oo 

Z= ^ n\n)(n\ +^n(|2n- l)(2n- 1| + |2n)(2n|) (18a) 

n=— oo n=l 

oo 

X+ = Y \n + l)(n\ + J2\n + 2)(n\ (18b) 

n=— oo n=l 

This configuration is not a static solution of the field equations for gauge theory on the fuzzy 
cylinder with action given by Eq. (12). Nevertheless, it demonstrates non-trivial topological 
properties one may expect for the domain wall solution. In NC field theory the role of 
topological invariants is played by traces of commutators. Indeed, the latter do not change 
under small variations of operators they depend on. For configuration (18) we have 

Tr[X + ,X_] = -l. (19) 

The existence of such topological charge suggests a strategy of the search for the domain 
wall solution: one may try to obtain a BPS bound for the energy functional following from 
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Eq. (12) in the topological sector determined by Eq. (19). However, as we will see below, the 
action (12) as it stands does not admit BPS solutions with desired properties — the domain 
wall tends to dissolve. 

To stabilize the domain wall, let us introduce the following additional term into the gauge 
theory action 

2 IPi T r (- — 
g 2 V 16/ 2 

This term has simple physical effect: it stabilizes the radius of the cylinder. In the commu- 
tative limit it becomes a mass term for the adjoint scalar field A p , 



S m = I dt -^Tr --—(X + X_ + X_X + - 2) 2 ) . (20) 



J dtpdQdz (^-^-A^j as I -> . 



With this term added, the static energy takes the following form, 

- 277 Tr (\ [Z, X + ] -X + \ 2 + £[X + , X_f + ^(X + X_ + X_X + - 2f ) , 21 ) 



g 2 pl 

where we introduced dimensionless parameters 



A=f. ,= ^. (22) 



This expression can be rewritten in the BPS form 
2tt 



g 2 pl 



Tr (\[Z,X + \ - X + | 2 + ± | ({X + ,X_} - 2)J 



=F^{[X + ,X_],({X + ,X_}-2)}) >T^^Q, (23) 
where braces stand for anti-commutators, and 

Q = ±Tr{[X + ,X-], ({X + ,X_} - 2)} = Tr[X + , (X_X + X_ - 2X_)] (24) 

is a topological charge. For the configuration (18) the topological charge is equal to one. 
Energy in each topological sector is minimized by the solution of BPS equations, which follow 
from Eq. (23), 

[Z, X+] -X + = (25a) 
[Z, X-] + X_ = (25b) 
(A + /i)X_X + ± (A - fi)X+X- = 2p . (25c) 
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These equations can be considered as defining an algebra with three generators Z, X + , X_. 
Each solution of Eqs. (25) decomposes into direct sum of operators acting in irreducible 
representations of this algebra. Thus, we have to classify irreducible representations of the 
algebra (25). 

Let us work in the eigenbasis {\z)} of the operator Z, 

Z\z) = z\z) 

It follows from Eqs. (25a), (25b), that X + , X_ raise and lower the eigenvalue z by one, 

X+\z) = x + {z)\z + 1) , X_\z) = x-(z)\z - 1) . 

Hence, in an irreducible representation, the eigenvectors of Z may be labeled by an integer 
k, and in this basis one has 

Z = J2(k + l)\k)(k\ 
k 

X + = J2 x k \k+l)(k\ 

k 

x_ = X>i !*><* + 1| , 

k 

where 7 is a real parameter characteristic to the irreducible representation. Relation (25c) 
gives a recursion equation for the coefficients x k . Let us consider the lower sign in Eq. (25c), 
then 

(A + n)\x k \ 2 - (A - n)\x k _ x \ 2 = 2/x . 
The solution of this equation is 

\x k \ 2 = 1 + Ca k 

where 

A — [i 
A + /I 

and C is an arbitrary constant. Taking into account that |a| < 1 and that \x k \ 2 > we 
obtain three different possibilities (up to redefinitions of k, C and a): 

a) C = and k runs from —00 to +00; in this case x k = 1, 

b) C = — 1 and k runs from 1 to +00; in this case x k = y/l — a k , 

c) C > and a > and k runs from —00 to +00. 

One may give the following geometric interpretation to these three types of solutions. 
Let us consider them as axially symmetric fuzzy manifolds where the role of coordinates is 



13 



played by covariant coordinates X + , X_ and Z. This kind of interpretation is quite common 
in NC gauge theories; later it will be supported by the study of the spectrum of small 
perturbations. The case (a) corresponds just to the original fuzzy cylinder (see Fig. la)). 
The case (b) corresponds to the semi-infinite fuzzy cylinder ("fuzzy fingerstall") of variable 
radius, starting at z = 7 and oriented towards positive values of z (see Fig. lb)). Finally, the 
case (c) corresponds to an infinite tube whose radius tends to unity as z — > +00 and becomes 
infinite as z — > —00 (see Fig. lc)). It is straightforward to check that the topological charge 





z 

a 




Figure 1: Geometrical interpretation of the three types of solutions to BPS equations (25) 
a) fuzzy cylinder, b) fuzzy fingerstall, c) infinitely expanding tube. 

Q is equal to zero in the first case, equal to one in the case (b) and is ill-defined in the case 
(c). We will not consider solutions of the third type in this paper. 

The solution corresponding to a domain wall between U(l) and U(2) gauge theories is a 
direct sum of representations (a) and (b), 







m\ 



Z = n\n)(n\ + J^n|2n- l)(2n- 1| + J^(n + 7)|2n)(2 

n=— oo n=l n=l 

00 00 

X + = ^ \n + 1) (n| + J2\2n + l)(2n-l\ + J2Vl- a n \2n + 2) (2 



71 



(26a) 
(26b) 



71=1 



n=l 



This domain wall has a simple geometrical interpretation: it describes a parallel system of an 
infinite fuzzy cylinder and a half-infinite fuzzy fingerstall (see Fig. 2). The topological charge 
Q of this solution is equal to unity, as expected. Without loss of generality we dropped in 
Eq. (26a) the possibility of the overall shift of Z by a c-number. Parameter 7, as we will 
show later, characterizes the position of the wall along ^-coordinate. The energy of the wall 
is 

nmp 



£ 



2g 2 P 



(27) 
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Figure 2: Geometrical interpretation of the domain wall: fingerstall inserted into the fuzzy 
cylinder. 



Its width 5z can be estimated as the size of the region along z-direction, where exact solution 
(26) differs considerably from its asymptotic form (17). Then, restoring the dimensionality, 
one has 

5z = r « — at m « - . (28) 

| In a | m I 

Note, that as m tends to zero, the width of the wall becomes infinite. This justifies our 
previous claim that in NC pure Yang-Mills theory, the wall tends to dissolve. Thus, we keep 
m 7^ 0. In the commutative limit I — > the energy of the wall diverges, while its width 
stays finite. Another interesting limit is that of the NC plane. It corresponds to p — > oo, 
I — > 0, and -& — pi fixed. In this regime the energy (27) diverges as p 3 , which means that our 
solution does not correspond to a domain wall of finite tension on the NC plane. 

The solution (26) is easy to generalize to (U(Mi) — [/(.A/2)) domain wall with Mi < A/2. 
One simply takes the direct sum of Mi representations of type (a) (cylinders) and M2 — Mi 
representations of type (b) (fingerstalls). The antiwall (that is (U(Mi) — U(M 2 )) wall with 
Mi > A/2) is the direct sum of irreducible representations of the algebra defined by relations 
(25) with upper sign in (25c). A system of a wall and an antiwall will be studied in detail 
in section 5. 



4 Scalars and fermions in the domain wall background 
4.1 Adjoint scalar 

To figure out the properties of the domain wall solution constructed in the previous section, 
let us consider a Hermitian adjoint scalar field 4> in the ([/(!)—[/ (2)) domain wall background. 
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The action is 



S^ = 2npl J dmQ(D o 0) 2 + ^([Z,0] 2 + [X,0] 2 + [Y,0] 2 )-im^ . 
After redefinition (16) one obtains 

S, = j J ^Tr(p 2 (D o 0) 2 + [^0] 2 -A 2 |[X + ,0]| 2 -^ 2 ) . (29) 
where = m^p. Equations for normal modes following from Eq. (29) are 

pV0 = [Z, [Z, 0]] + y[X_, [X + , 0]] + y [X + , [X_, 0]] + // 2 , (30) 

where operators Z, X + are given by (26). It is convenient to decompose the Hilbert space 
H into direct sum H = H 1 © H 2 of two Hilbert spaces with bases 



\c n ) 



\n) n < 

\2n-l) n>\ 



and 

\h p ) = \2p), P >1 
respectively. Background operators (26) take the following form 



Z= ^ n\c n ) (c n \ + ^(p + -f)\h p ) (h p \ (31a) 

n=— oo p=l 



oo oo 



X + = J2 \cn + i){c n \ + ^^l-aP\h p+1 )(h p \ . (31b) 

n=— oo p=l 

These operators are block-diagonal in the sense that they map Hi to Hi (i = 1,2). It is 
convenient to decompose the field as 

= ^2<j>nm\c n }(Cm\ + ^ tfipqlK) (K\ + ^(Xnp\c n ) (hp\ + Xnp\ h P ) ( C «l) » ( 32 ) 

so that equations for different components 0„ m , and x„ p decouple. One immediately 
notices that the indices p, q are greater than zero, so nonvanishing matrix elements of the 
operators ip and \ have at least one positive index. This implies that the modes described 
by these operators live in the region z > and cannot penetrate into the region z — > — oo. 
Let us study the three types of modes entering the decomposition (32) separately. 
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Equations for the nm -components in the background (31) exactly coincide with equations 
for the scalar field on the fuzzy cylinder (9). Thus, excitations of this type do not feel the 
presence of the fingerstall and freely propagate from one end of the cylinder to the other. 

Equations for the (^-components have the following form, 

22 (, x2 2 l2 f a p + a p ^ + a q + a q - 1 \\ 
pWip pq = (p - qf + /A + 2A 2 1 <p pq 



4 

- AV(l-a p )(l-a«) " AV(1 - - a 9 " 1 ) ¥> P -i,g-i • 

Let us perform Kaluza-Klein decomposition in analogy to the case of free scalar field con- 
sidered in section 2. Namely, after fixing the number of the Kaluza-Klein mode, p — q = N, 
we obtain the following recursion relation for ip q = <p q+Njq , 



p a; y? 9 = I N + /i + 2A II 11^ 

- AV(1 " a N +i)(l - a q ) Vq+1 - AV (1 - e^+^Xl - a'" 1 ) ^-i • (33) 
At large q, the general solution of this equation is 

Vq = Ae lklq + Be~ iklq , (34) 

where the wave vector k and frequency u are related by Eq. (8). To relate the coefficients A 
and B let us note that equation (33) leads to the conservation of "current", 

Jq+l — Jq = 

with 

J q = i^(l-a N +«)(l-a q ){^ q+lVq - tp* q <P q+ i) ■ 
Equation (33) with q — 1 implies that 

with real coefficient of proprtionality c. So J\ = and, as a consequence, we obtain that 
J q = for all g. Then the coefficients in Eq. (34) have equal absolute values, 

\A\ = \B\ . 

The asymptotic solution (34) describes two waves with equal amplitudes propagating in 
the opposite directions. In other words, excitations f pq experience total reflection from the 
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domain wall. Note that these modes do not feel the presence of the fuzzy cylinder and live 
entirely on the fingerstall. 

Modes Xnp obey the following equation, 

p 2 ^ 2 Xn P =({n-p- 7 ) 2 + Hi + 2A 2 (\ - aP + ^ P ^ ^ Xnp 

- AVl - aP- l Xn-i, P -i ~ AVl - « p Xn+i,p+i • 
By fixing n — p = N and substituting x P = X p +n, p one obtains 

pV Xp = (V - 7 ) 2 + 2A 2 (l - ^^)) X P 

- AVl - aP- 1 Xp-i " AVl - • 

The analysis similar to that for modes ip demonstrates that solutions of this equation are 
waves which experience total reflection from the domain wall. Their dispersion relation is 

2 (iV- 7 ) 2 (2 kl\ 2 2 

^= y —-2L+(-sm-) + m 2 (35) 



P 2 V 2 

Let us comment on this result. First, if the parameter 7 is not integer, spectrum of the 
off-diagonal modes Xn P differs from that of diagonal modes (p nm and ip pg . As z — > +00, 
the x an d (0, V 9 ) modes correspond to off-diagonal and diagonal components of the adjoint 
U(2) field, respectively. Difference in their dispersion relations means that in the case of 
non-integer 7 the gauge group U(2) is broken down to U(l) x U(l) at z — > +00. This effect 
is nothing else than spontaneous breaking of the gauge symmetry by a nontrivial Wilson line 
on the cylinder. Indeed, from (13c) and (26a) we see that 

A e = (8°) , . 

The value of the Wilson line around the cylinder is equal to e 27 ™ 7 7^ 1 if 7 is not integer. 
Thus, Ae cannot be removed by a gauge transformation in this case and breaks the gauge 
symmetry. Conversely, if 7 is integer, Aq is a pure gauge as z — > +00, and the symmetry 
[/ (2) is unbroken in the z — > +00 asymptotics. From now on we consider integer values of 7 
only. 

Let us see now that the parameter 7 is related to the position of the domain wall. Equa- 
tion (35) implies that the genuine Kaluza-Klein number of a x excitation which characterizes 

its energy is N = N — 7. The expression for iVoth KK excitation in the operator form is 

00 

X (N 0) ^J24 N0) \C 1+ N 0+P )(h p \ 
p=l 
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The operator x^ N °Hx^ N °^Y ac ts diagonally in Hi and annihilates all vectors with 

n < n min = N + 7 + 1 . 

Hence, the x _wave is reflected from a point with ^-coordinate n min . For a given KK mode 
this coordinate depends additively on 7, implying that this parameter is natural to interpret 
as the position of the wall 5 . 

Let us summarize our results on the spectrum of the adjoint field in the presence of the 
domain wall. In the region z — > +00 there are four (real) modes which constitute adjoint 
representation of U(2). Their dispersion relation is given by Eq. (8) which reduces to the 
standard dispersion relation on the cylinder in the long wavelength regime. Three of these 
modes are reflected from the wall while one freely propagates to the region z — > — 00. A 
notable fact is that the latter mode does not decouple from the other modes in the region 
z — > +00: in terms of U(2) group it has the form q), and couples to other modes due 
to gauge interactions. These results can be straightforwardly generalized to the case of 
the (U(Afi) — U(J\f 2 )) domain wall with A/i < A/" 2 . In that case there are A/"f modes in 
the adjoint of U{Mi) group freely propagating along the cylinder and (A/" 2 2 — A/" 2 ) modes 
exhibiting reflection from the domain wall. 

The analysis of gauge field perturbations in the domain wall background is less trans- 
parent because of mixing between different components and the necessity to implement the 
Gauss' constraint. We do not give the details here, and only states that the qualitative 
picture described in the previous paragraph holds for the gauge field perturbations as well. 
Namely, the U{Mi) gauge fields propagate freely along the entire cylinder; in the region 
z — > +00 they become a part of £7 (A/2) gauge multiplet; the rest of the U(Af 2 ) gauge fields 
live only in the region on the right and are totally reflected from the wall. 

4.2 Fundamental fermion 

To study fermions in the background of the domain wall, we have to define the Dirac operator 
on the fuzzy cylinder. To this end let us consider possible choices of the Dirac Hamiltonian 
on the commutative cylinder. The simplest one is 

D = ia 2 d z a 3 do + <j\mt 

P 

5 Strictly speaking, the above argument is not rigorous. The reason is that bilinear combinations of adjoint 

fields are again adjoints and are not gauge invariant, while in NC theory gauge transformations of adjoints 

include change of coordinates. Our conclusion that 7 characterizes the position of the domain wall will be 

confirmed by the consideration of a fundamental field in the next subsection. 
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where to/ is the fermion mass and Oi are Pauli matrices. This operator is not convenient 
for our purposes, because natural derivative operators on the fuzzy cylinder are xd z and yd z 
rather than d z (see Eqs. (5a), (5b)) 6 . So, we will work with the following unitary equivalent 
operator 

D = S+DS , 



where 

' iS/2 Q 



(36) 



e- 10 / 2 
Explicitly, one has 

D = ^ (^-i{<Jiy - o 2 x)d z - ia 3 dg + ^ + (a ± x + <j 2 y)m^j . (37) 

Let us note that due to the form (36) of the unitary transformation S, periodic boundary 
conditions for the operator D correspond to anti-periodic boundary conditions for D and 
vice versa. This is not problematic, as one may consider both periodic and anti-periodic 
spinors on the cylinder. 

Noncommutative counterpart of the operator (37) can be naturally defined in terms of 
its action on a spinor field tp as (cf. Eqs. (4), (5)) 

D^j = -(<7i[x, V] + <r 2 \y, if>] + <r 3 [z, ^\) + — V> + ^(<ri{x, ^} + a 2 {y, ^}) . (38) 
pi 2p 2p 

The operator ordering in the mass term in Eq. (38) is chosen somewhat arbitrarily; a different 
choice would not alter the results of the analysis below but would make calculations more 
cumbersome. The spectrum of the operator (38) is 

2 (iV + 1/2) 2 (1 . kl\ 2 / A;/\ 2 



p 2 1 yi 2 ; 1 v 1 2 



6 Nevertheless, one could proceed with the operator D and rewrite the term ioid^ as 

i-^(x ■ xd z 4> + y ■ yd z ip) 
P 

Noncommutative counterpart of this expression would be 

-^j (x[y , ip] - y [x, V>] ) = {xtpy - y^x) 

This operator leads to the second order recursion equations and its spectrum suffers from doubling of fermion 
species, resembling that occurring in the lattice field theory. 
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where N is the Kaluza-Klein number and k is the wave vector. Two comments on this 
dispersion relation are in order. First, the appearance of (N + 1/2) instead of N in the first 
term of Eq. (39) is not unexpected, since as discussed above we are effectively considering 
anti-periodic fermions. The dependence on k in the third (mass) term is more surprising. 
Interestingly, the same dependence appears in the dispersion relation for the field A p defined 
by Eqs. (13a), (13b) when the term (20) is added to the gauge action. 

The Dirac Hamiltonian for fermion in the fundamental representation of the gauge group 
is obtained from Eq. (38) via substitution 

[x, if)] — > — ^x 
{x, ijj} — > X"0 + 

and similarly for other coordinates. In the domain wall background (26), it is convenient to 
decompose the fermion field in the following way (cf. Eq. (32)) 

oo oo 

^= ^2 ^nm\c n )(m\+ ^ Vpn\h p )(n\. 

n,m=— oo p=l,n=— oo 

Equations for the ?/> nm -components are identical to those obtained in the case of fuzzy cylinder 
without the domain wall, and the corresponding branch of the spectrum is given by Eq. (39). 
Thus, these modes do not feel the presence of the domain wall at all. For spinors r\ vn we 
obtain the following equations, 




(p-n + 7 + \W pn + (A + ^)VT^V 2 P+ i,n - (A - f )<„_! 
(X + ^VT^^vUn - (A - f - (p - n + 7 - \) ni 



2 

pn / 



where fj,f = mfp. After substituting rj^ = r}p~, p +Ni = V%, p +N-n equations along diagonals 
n — p = N take the form 

(-N + 7 + \H + (A + f )VT^V 2 P+1 - (A - f W P 




+ f )VT^ rT V 1 P -i - (A - f W P ~ (-N + 7 + \K, 

Since n = p + N > N, a KK mode with fixed iV does not penetrate into the region z — > — oo 
and thus is reflected from the domain wall. In this respect the situation is similar to the 
case of the adjoint scalar considered in the previous subsection. 

Let us find now a zero mode localized on the domain wall. By zero mode we understand a 
mode which is annihilated by the transverse part of the Dirac Hamiltonian. In other words, 
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it obeys the following equations 

(A + f ) VT^^vU - (A - f ) vl = (40) 
(a + f ) VT^v 2 P+1 - (a - f ) ^ = . (41) 

Equation (40) with p = 1 implies that 77J = 0, and as a consequence 

^ = (42) 
for all p. From Eq. (41) one obtains the following solution for 77^-components 



p-i 

vl^C^Hil-a^ , (43) 
3=1 

where 

X-jii/2 
P X + fif/2' 

and C is a normalization constant. This solution is well-behaved at large p for jif > 0. Note 
that zero mode (42), (43) is chiral 

vzV = -V ■ 

A zero mode with opposite chirality is localized on the antiwall configuration. 

Let us work out the profile of zero mode along z-direction. This profile is gauge invariant 
and can be regardedd as a probe of the domain wall shape. We assume fx, \x$ A. First, let 
us note that the energy of zero mode on the iV-th diagonal is given by 

N-j-1/2 
uj = . 

P 

Thus, the genuine Kaluza-Klein number of this mode is Nq — N — 7 — 1. From Eq. (43) one 
obtains the following gauge invariant density 

00 

rf r] = J2(ri 2 p ) 2 \N + -f + p)(N + 1 + p\. (44) 
P =i 

Its Weyl symbol is shown in Fig. 3. This density concentrates near p = p ma x with p max 
determined by 

1 - a Pmax = P 2 

or, explicitly, 

v ™ = T^Ju7 f (45) 
22 



z max 



Figure 3: Gauge invariant profile of zero fermion mode localized on the domain wall. 

Coefficients entering Eq. (44) may be approximated in the vicinity of p ma x as follows, 

(vlT = & exp (-^(P - Pmax?) (46) 

where C is a constant. From Eqs. (44) and (46), one infers the following properties of the 
Weyl symbol rfr)(z) in terms of physical parameters: it is equal to zero at z < z min , where 



peaks in the vicinity of 



where it has the form 



h + IN , (47) 



-ln_L- + Z 7 + Z7Vo. (4 8 ) 
to 2m f I 



rfr)(z) = C 2 exp (-m f m(z - z max ) 2 ^j , (49) 

and falls off exponentially at z — > +oo. The coordinates z min and z max can serve as two 
alternative definitions of the position of the domain wall. If N is not too large, the last 
terms in Eqs. (47), (48) can be neglected. Then, z min ~ Z7 coincides with our previous naive 
estimate of the domain wall position, based on gauge dependent quantities (see discussion 
after Eq. (35)). 

On the other hand, z max is more appropriate for the definition of the position of the 
wall, if one considers rfrj(z) as the gauge independent shape of the domain wall seen by the 
zero fermionic mode. The disadvantage of this definition is that it depends on the mass of 
the fermion, but it is problematic to provide a probe independent meaning to the notion of 
the domain wall shape. Assuming to/ ~ to, one obtains from Eqs. (48), (49) the following 
estimates for the position z max and width 5z of the domain wall 

Zmax ~ — In —^7 + h , 5zk, — . (50) 
to 2ml m 
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Let us stress once more that these estimates refer to the domain wall as seen by the zero 
fermionic mode, and other probes may, in principle, give other results. The estimate for the 
width agrees with that deduced from the (gauge dependent) profile of the domain wall, see 
Eq. (28). 

5 Wall-ant i wall system 

Similarly to the domain wall solution, it is convenient to describe the wall-antiwall configu- 
ration in terms of operators acting on a direct sum of Hilbert spaces H = HiQ)H 2 . Subspaces 
Hi and H 2 emerging in the case of the wall-antiwall system are spanned by the following 
systems of basis vectors 

\c n ) = \n) , n < 

\K) = b> » v > i ■ 

In these notations the wall-antiwall field has the following form 

oo 

z = ^2 ( n + 7i)|c n )(c n | + ^2(p + l2)\h p )(h p \ (51a) 

n=— oo p=l 

oo 

X + = Vl - a- n \c n+1 )(c n \ + - aP \ h P+i)( h p\ ( 51b ) 

n=— oo p=l 

This may be viewed as the union of two fingerstalls. Two physically different situations 
occur depending on the value of the parameter A7 = 72 — 71: for A7 > the two fingerstalls 
intersect (see Fig. 4a) and, as we will see below, U(2) gauge theory emerges in the region 
between the walls, while for A7 < the fingerstalls are disconnected (Fig. 4b). 

Due to the direct sum structure of Eqs. (51), they describe an exact solution of field 
equations. This is somewhat unexpected, because it implies that there is no attraction 
between the wall and antiwall. On the other hand, configuration (51) belongs to the trivial 
topological sector (it is straightforward to check that its topological charge Q defined in 
(24) is equal to zero), and, consequently, it is expected to be unstable. We will show this 
explicitly later on. 

To justify our interpretation of the two-fingerstall configuration shown in Fig. 4a as a 
wall-antiwall system, let us address a question of the restoration of non-Abelian symmetry 
in the region where the fingerstalls overlap. When the size of this region is large, 

A 7 » 7T~~ r , (52) 
men 
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Figure 4: Two possibilities for the non-BPS solution of two fingerstalls: a) Two regions with 
U(l) theories separated by a region with U{2) theory; b) Two disjoint regions with £7(1) 
theories 

we expect U(2) gauge symmetry to emerge in the overlap region. Direct analysis of gauge 
field fluctuations in the background (51) is rather cumbersome. More illuminating, as in the 
case of the single domain wall, is to consider an adjoint scalar field 0. It obeys Eq. (30) 
with the background given by Eqs. (51). For simplicity we consider massless field, = 0. 
In analogy to the single domain wall case, we perform the decomposition (32). Dispersion 
relations for modes <f) nm and ip pq are again given by Eq. (8). These are modes living on each 
of the fingerstalls separately. Modes <f) nm do not feel the presence of the domain wall, but are 
reflected back from the antiwall, while modes <p pq do not feel the antiwall but are reflected by 
the wall. The same qualitative picture is valid for the modes of the gauge fields responsible 
for U(l) gauge symmetry in the asymptotic regions z — > — oo and z — > +oo. This means 
that the gauge fields from z — > — oo region cannot penetrate into z — > +oo region, and vice 
versa. 

Modes <p nm and ip pq correspond to diagonal elements of the U (2) multiplet in the region 
between the walls. Off-diagonal components of the multiplet should come from the Xnp 
sector. Equation for these modes reads 

p 2 ^ 2 Xn P = l(p-n-A 1 ) 2 + 2\ 2 (l J J Xnp 



- A 2 a/(1 - aP-!)(l - a-" +1 ) Xn-iv-i ~ AV(l-« p )(l-a-")Xn + i, P +i • 

Performing the Kaluza-Klein decomposition, p — n = N, Xp+n, p = Xp > one obtains 

92 ( , t v, ,9 { aP + aP" 1 + a N ~ p + a N - p+1 \\ 
pW Xp = (iV-A 7 ) 2 + 2A 2 1 ))xp 



- \ 2 y/(l- a p-i)(l- a N-P+i) Xp -i ~ \ 2 ^(l-aP){N-p) Xp+ i ■ (53) 
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The index p in Eq. (53) runs from 1 to AT because n < 0. If the condition (52) is satisfied, 
low energy modes are those with 

|JV | = \N - A7I < A7 

and their dispersion relation is approximately given by Eq. (8) (with N instead of N). 
Strictly speaking, the wave vector k of these modes is not arbitrary and takes discrete values 
determined by the boundary conditions at p — 1,N. But the gap between two neighboring 
values is of order 

Ak 



/A 7 

and tends to zero as A7 — > 00. Thus, in the limit of large separation between the walls, the 
spectra of the off-diagonal and diagonal modes coincide, and U (2) gauge symmetry emerges 
in the region between the walls. 

Now, let us return to the question of stability of the wall-antiwall solution. As there is 
no force between the wall and antiwall, instability, if any, must reveal itself in the existence 
of a tachyonic mode in the spectrum of fluctuations of gauge fields about the wall-antiwall 
background. Clearly, this tachyon should carry indices corresponding to both fingerstalls 
comprising the solution (51). Thus, we consider fluctuations of the form 

Z = Z {0) - a e , X + = xj 0) + b + 

where Z^°\ are given by Eqs. (51), and a e , b + are off-diagonal 

clq = ^^(a pn \h p ) (c n \ + a* n \c n ) (h p \) 
b+ = 'y^ J (b pn \hp) (c n \ + c np \c n ) (h p \) . 

Insertion of these expressions into Eq. (21) yields quadratic energy functional for the fluc- 
tuations. The latter can be expressed as the sum of contributions coming from different 
diagonals p — n — 1 = N. The contribution due to the A^tli diagonal is (we drop the overall 
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dimensional prefactor in the energy) 



p=N+l 



£n= [ ( N - A ^ b P - V 7 ! - « iV -P+ 1 a p + Vl - a^ 1 a p _i 



p=i 



+ (AT — A7)c* + Vl - aP- 1 a p - \/l - a N -P+ l a p _i 



A 2 

4 



Vl - a p - x c; + v 7 ! - a w " p+1 & P - \fl-a N -Pc* p+1 - y/l - aP b 



'p+i 



Vl - a^ 1 c* + v 7 ! - a w " p+1 & P + V 7 ! - ct N ~P c* p+1 + v 7 ! - « p 6. 



'p+i 



+ ^(a^- 1 + a N ~ p+1 )\c p \ 2 - ^(aP^ 1 + a N - p+1 )\b p \ 



(54) 



Notations in this expression are 



b p ,-N+p-i 



a, p — a,p-N+p 



Cp C_jv+p,p— 1 



p = l,2, ... ,N+1 
p = l,2, ... ,iV 
p = 2,3, ... ,iV. 



The expression (54) is rather lengthy, and we presented it just to demonstrate that the only 
negative contribution to the energy comes from the modes b p . This suggests that these modes 
give the dominant contribution to the tachyon. 

The following analysis is different for the two cases shown in Fig. 4. Let us consider first 
the case A7 < 0. We argue that in this case the lowest eigenvalue of the energy is given 
by the contribution of the zeroth diagonal. Indeed, modes along other diagonals have larger 
Kaluza-Klein energy (which is proportional to (N — A7) 2 and N > 0). This effect shifts 
upwards masses squared of these modes. Thus, let us study the case N = 0. The expression 
(54) is greatly simplified yielding 



S = ((A7) 2 - 2/j,\ S )\b 1 \ 2 . 



We see that 61 is an eigenmode with the frequency 

(A7) 2 



m 
T 



This means that at A7 < 0, the system is unstable, provided that IA7I < p\/m/l, but 
becomes stable at large values of IA7I. This behavior is fairly natural, because when A7 < 
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and IA7I is large, two fingerstalls comprising the wall-antiwall system do not intersect and 
are well separated (Fig. 4b). 

Let us now turn to the case A7 > 0. If A7 satisfies Eq. (52), there is a large region of 
intersection of the two fingerstalls with U(2) gauge theory inside (Fig. 4a). One expects the 
system to be able to roll down to either U(l) or U (2) vacuum on the whole cylinder. Instead 
of trying to find the tachyonic mode exactly, let us present an Ansatz which demonstrates 
that the energy functional given by Eq. (54) has negative directions. To this end, let us 
choose N = A7, take a p = 0, c p = and rewrite the energy £at=a 7 in the following form 

N/2 

£n=Aj = 



+ 



2 ( ( A + ^ a Vl-a N - p+1 b P - VT^7^b p+1 - A\^ia N - p+l \b p \ 2 
p=i ^ 

N+l . 



p=7V/2+l 



V 7 ! - a N ~P+ 1 b p - aVl-aPh 



AXfxa p 1 \b p \ 



(55) 



-2X^(1 - a N / 2 )\b N/2+1 \ 
Positive terms in eq. (55) can be set to zero by the following choice 



.7=1 v 



b P = < 



a 



p-i 

N+l-p J-J 
3=1 



p=l, ... ,N/2 



p= N/2 + 1, ... ,N + 1. 



(56) 



This Ansatz for b p is a symmetric combination of two bell-shaped functions localized around 
Pmax and (iV — Pmax), where p max is given by Eq. (45) (with \x instead of ///). Substitution of 
the Ansatz (56) into the energy functional (55) yields the following estimate for the tachyon 
energy 



u 2 = -C 



■ a 



A 7 -2 Pl , 



= —Cm e 



-mAz 



where Az = /(A7 — 2p max ) is the separation between the wall and antiwall, and C is a 
coefficient of order one. We see that the tachyon mass is exponentially small when the walls 
are far away from each other. 

This tachyon is not directly related to the distance between domain walls A7 which is an 
exact modulus of the wall-antiwall solution. One may guess that tachyon condensation leads 
to the change of the shapes of the walls. The natural candidates for the end points of the 
tachyon condensation are U(l) and/or U(2) vacua. We leave the study of this condensation 
for future. 
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6 M(atrix) theory interpretation of the domain walls 



The purpose of this section is to suggest a way of embedding the domain wall solutions 
constructed in this paper into the matrix model of M-theory (see, e.g., Ref. [12] for a review 
of the M(atrix) model). Matrix model is supersymmetric quantum mechanics described by 
the following Lagrangian 

L = ^Tr jx* + ipr.X'"] 2 + (fermions) j , (57) 

where X* (i — 1, . . . , 9) are real-valued N x N matrices subject to the constraint 

[X 1 , X*] = . 

Originally [13], this quantum mechanical system was suggested as a regularized theory of 
a (super)membrane in flat 11-dimensional space-time in light-cone gauge. The matrices X* 
play the role of embedding functions of the membrane. Regularization is removed by taking 
the limit N — > oo. In this language R = 27r/f 1 is the membrane tension. 

Alternatively, one may consider Lagrangian (57) as an effective low-energy description 
for a system of N DO-branes in the type IIA-theory in the A = gauge. In this case R has 
an interpretation of compactification radius of 11-dimensional M-theory to ten dimensions 
and in string units, l s — 1, this radius is equal to the string coupling g s . 

It was conjectured in Ref. [14], that large- N limit of the matrix model (57) describes all 
of the M-theory in the infinite momentum frame. Moreover, it was argued [15, 16] that the 
quantum mechanical system (57) at finite N describes a sector of M-theory with N units of 
momentum along compact light-like direction. 

In Ref. [8] fuzzy cylinder was obtained as a BPS-solution in the matrix model (in A = Z 
gauge) and was interpreted as a D2-brane of type IIA theory. In Ref. [11] field configurations 
similar to our Eqs. (18) were discussed in this context and were interpreted as junctions 
of D2-branes. It is worth noting, however, that these junctions (domain walls) were not 
obtained as solutions of matrix model equations. 

Here we would like to suggest that domain walls studied in this paper may be obtained 
as solutions of the matrix model in curved backgrounds. 

If one sets X* = for i > 4 in the matrix model Lagrangian, one arrives at the action 
very similar to the action (12) of gauge theory on the fuzzy cylinder. The only difference 
is that extra terms linear in X and Y present in the definition of the field strength in this 
theory (see Eqs. (lib), (11c)) are absent in the matrix model. So our purpose in this section 
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is to find a way to introduce these terms, as well as the term (20), into the matrix model 
Lagrangian. 

A generalization of the matrix model Lagrangian (57) to arbitrary curved background is 
not known (see Ref. [17] for a discussion of this problem). However, there is a proposal [18] 
on how to modify the Lagrangian of the matrix model to incorporate the effect of arbitrary 
weakly curved background independent of the light cone coordinate x~ . Namely, to describe 
the effect of non-trivial eleven-dimensional metric gMN = Vmn + h>MN and three-form Amnl 
at the linear level, one adds the following terms to the Lagrangian of the matrix model 

AL f = xE^--^STV(T MW X il ...X'") (58) 
AL A = ^2-d il ...d in A MNL STr(J MNL X il ...X i -) (59) 

71/. 

where STr stands for the totally symmetrized trace 



STr(A ai ...A an ) = - ] TiJ2 A a{ai) ...A 



transmutations <7 

We have not written terms describing magnetic interactions of the membrane and terms with 
fermions. The former do not appear in the backgrounds considered below, while the latter 
are not relevant for our purposes. 

Explicit expressions for the components of the energy-momentum tensor T MN and anti- 
symmetric current J MNL can be found in Ref. [18]. In what follows we will make use of the 
expressions for the T ++ and J +tj components 

T ++ = i (60) 

J +iJ = ~\x.\x']- 

Eqs. (60) and (58) suggest that a natural starting point to construct the M-theory back- 
ground leading to extra terms like (20) in the matrix model action, is to consider metric 
with non-trivial g ++ component. In order to be a legitimate background of the M-theory at 
least in the supergravity approximation, this metric should be supplemented with the ap- 
propriate three-form field to satisfy the equations of eleven-dimensional supergravity. Thus 
one naturally arrives at the following class of the supergravity solutions [19] 

ds 2 = -2dx + dx~ +^dx l dx l - H(x l )(dx + ) 2 

F+ijk £,ijki% ) i 
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where functions H{x l ) and (,ijk(x l ) are related as follows 

%h = l -i ijk e jk (6i) 

and 

£ = -iijkdx 1 A dx j A dx k 

is closed and co-closed form. These solutions are generalizations of the homogeneous pp-wave 
solutions [20] which have attracted much attention recently. 

To start with, let us present an example of the background admitting fuzzy cylinder as a 
solution of field equations. For this purpose it suffices to consider quadratic function H(x l ). 
Say one considers H(x l ) of the following form 

H fc (x i ) = l 2 ((x 1 ) 2 + (x 2 ) 2 ) (62a) 

supplemented with the following three-form 

£ fc = lduj , u = (x l dx 2 A dx 3 - x 2 dx 1 A dx 3 ) . (62b) 

Applying the rules described above, it is straightforward to check that the bosonic part of 
the matrix model Lagrangian in this background has the following form 

L fc = ^-Tr (X*X< + ([X 1 , X 3 ] + tlX 2 ) 2 + ([X 2 , X 3 ] - ilX. 1 ) 2 + [X 1 , X 2 ] 2 + . . . ) , (63) 
2R V / 

where dots stand for non-negative terms in the potential, vanishing for zero X 4 , ...,X 9 . 
One immediately recognizes that the expression (63) is the same as the Lagrangian of the 
gauge theory on the fuzzy cylinder where X X ,X 2 ,X 3 play the role of covariant coordinates 
X,Y, Z and X 4 ,...,X 9 are massless adjoint scalar fields with a specific positive definite 
quartic potential. 

To construct gravitational background leading to the matrix model admitting domain 
wall solutions discussed above, let us note first that for two pp-wave solutions described 
by functions H^ 2 and three-forms £1,2, the sum H\ + H 2 and £1 + £2 is again a solution, 
if there are no hyperplanes with non-vanishing fluxes for both three-forms £1 and £ 2 - We 
have already described gravitational background leading to the gauge theory on the fuzzy 
cylinder, so now we have to find the function H m and three-form £( m ) such that £( m )i23 = 0, 
which give rise to extra term of the form (20) in the matrix model action. 

A natural guess would be to take 
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However, the Laplacian of this function is not positive definite, 

in contradiction to Eq. (61). 

To get around this difficulty, one may consider background depending on larger number 
of coordinates. For instance, one may make the following choice of the function H m and 
three-form £( m ) 

£( m ) = 2-\/2 y duj , uj — (x 1 x 4 dx 2 A dx 5 + x 2 x 5 dx 1 A dx 4 ) 
A 

The extra piece in the matrix model potential coming from this background is 

V m = |^STr ((X 1 ) 2 + (X 2 ) 2 - l) 2 + . . . (64) 

where dots stand for the terms which are at least second order in coordinates X 4 , X 5 . These 
terms do not affect equations for configurations with X 4 = X 5 = which we are focusing on 
here. Now, it is straightforward to check that 

^STr ((X 1 ) 2 + (X 2 ) 2 - I) 2 = ^Tr ((X 1 ) 2 + (X 2 ) 2 - l) 2 + ^Tr[X\ X 2 ] 2 + . . . (65) 

where dots now stand for commutator terms which do not affect the field equations (but, 
in the N — > 00 limit, in general contribute to the energy). The first term in Eq. (65) 
coincides with the extra term given by Eq. (20) while the second one can be eliminated by 
the redefinition of the parameters /j, and A (at /i 2 < A 2 ). Thus, Eq. (65) provides an example 
of the supergravity background leading to the matrix model with domain wall solutions 
discussed in this paper. It is clear from the discussion above that one can construct a variety 
of backgrounds with this property. 

There is a subtlety conserning the stability of the domain wall solution in the matrix 
model. As we showed in section 3, if the extra coordinates X 4 , ...,X 9 are disregarded, 
the domain wall saturates the BPS bound (23), and thus it is stable. However, when the 
extra coordinates are included, the omitted terms in Eq. (64) may become tachyonic, in the 
domain wall background, along the coordinates X 4 , . . . , X 9 , so that the stability may be lost. 
In order to ensure the stability, it is desirable to find a background admitting the domain 
wall solution, which saturates a BPS-type bound for the full matrix model action. 
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Such background can be obtained using technique proposed in Ref. [21]. The following 
class of pp-wave backgrounds with four extra 7 supersymmetries was described there, 



H 



dW(<j>) 



(66a) 
(66b) 



where a = 1, 2, 3 and 



<P 1 = x l + ix 4 
(f) 2 = x 2 + ix 5 



x 6 + ix 7 



Also, it was suggested in Ref. [21] that the matrix model potential for this background has 
the following form 8 



V* = Tr ( -[$ a ,$ 



+ 2 



^6 abc [$ b ,<S> c ]+d a W(<!>) 



(67) 



where <3> a , a = 1, 2, 3 are matrices corresponding to the coordinates <p a , and superpotenial 
W(&) is defined as follows, 

W($) = STrW(<Z>) . (68) 

For quadratic superpotentials the ordering prescription following from Eqs. (67), (68) agrees 
with that defined by Eqs. (58), (59) 9 . For generic superpotential, the two ordering pre- 
scriptions are different. We have nothing to say about this discrepancy here. If one adopts 
ordering defined by Eqs. (67), (68), background siutable for our purpose is a sum (in the 
sense explained above) of the pp-waves described by (62) and the pp-wave given by (66) 
with superpotential 

WrM = ty> ((0 1 ) 2 + (0 2 ) 2 - 1) • 



7 Generally, pp-wave background leaves unbroken 16 supersymmetries out of 32 supersymmetries present 
in the eleven-dimensional supcrgravity. These supersymmetries disappear after gauge fixing of the kappa- 
symmetry in the supermembrane action, so the resulting matrix model is not supersymmctric. In some cases 
pp-wave may have extra unbroken supersymmetries, and the corresponding matrix model is expected to be 
supersymmetric. 

We present terms in the potential which do not contain X ^ X •) X 9 coordinates. The latter terms arc the 
same as in flat space. 

9 In particular, by choosing W = l 2 ((0 1 ) 2 + {<p 2 ) 2 ) one may obtain supersymmetric realization of the fuzzy 
cylinder in the matrix model (with coordinates on the fuzzy cylinder X = Rc$i, Y = Re$ 2 , Z = Rc$ 3 
and all other fields set equal to zero), which apparently differs from that of Ref. [8]. 
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It is straightforward to check that in this background the domain wall (26) emerges as BPS 
solution (with identifications X = Re $ l5 Y = Re $2 an d Z = x 3 and all other coordinates 
set equal to zero). It is worth noting, however, that (2, 2) supersymmetry is explicitly broken 
in this background by terms coming from Eqs. (62). 

7 Summary and discussion 

Let us summarize our results and discuss some open problems. In this paper we constructed 
and studied domain walls between vacua with different gauge groups U{H\) and [/(.A/2) 
on one of the simplest NC manifolds, fuzzy cylinder. We demonstrated that these domain 
walls are characterized by a non-trivial topological charge and satisfy BPS-like equations, 
provided an extra term stabilizing the radius of the fuzzy cylinder is added to the action. 
They represent a novel class of exact NC gauge solitons in the sense that they cannot be 
obtained by making use of solution generation technique of Refs. [23, 24, 25, 26], and do not 
have commutative counterparts. 

By making use of the adjoint scalar and fundamental fermion fields as probes, we studied 
some of the properties of the domain walls and demonstrated that these objects exhibit rich 
pattern of non-trivial phenomena. Namely, we addressed a question, whether fields charged 
under the gauge groups can penetrate from one side of the wall to the other. The result 
is that if A/i < A/2, fields charged under [/(A/i) freely penetrate through the domain wall 
into the [/(A/2) region where they become a part of [/(A/^) multiplet. On the other hand, 
[/(A/2) fields which are not part of the [/(A/i) subalgebra experience total reflection from the 
domain wall. It is worth mentioning here, that the higher the mass (KK number) of the 
[/(A/2) mode, the deeper it penetrates into the region with the [/(A/i) vacuum. This effect 
is an illustration of the UV/IR mixing characteristic to NC theories. 

For fermion field we found that there is a zero mode localized on the domain wall. Wave- 
function profile of this mode is a gauge invariant characteristic of the shape of the domain 
wall. It would be interesting to understand whether there is a NC analogue of the index 
theorem relating the existence of this zero mode to non-trivial topological properties of the 
domain wall. 

Also we studied a non-BPS wall-antiwall configuration, dividing the cylinder into three 
regions with U(l), U(2) and again U(l) gauge theories. A somewhat unusual property of this 
system is the absence of the interaction potential between wall and antiwall at the classical 
level. Still, we found a tachyonic mode in the spectrum of perturbations about this system, 
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which we expect to roll down to either U(l) or U(2) vacuum on the whole cylinder. The 
precise mechanism of the tachyon condensation deserves further study. The mass of the 
tachyon becomes exponentially small at large separation between the wall and antiwall. 

There is a simple brane picture for the domain walls between £/(A/"i) and {/(A/^) theories. 
Namely, one can think of this configuration as a stack of A/2 — M\ branes of fingerstall shape 
inserted into Ai cylindrical branes. Some of the properties discussed above have a natural 
interpretation in the D-brane language. For instance, the pattern of gauge field spectrum in 
the presence of the domain wall is quite transparent in the string context. Indeed, gauge fields 
belonging to the £/(Ai) multiplet correspond to strings with both ends on the cylindrical 
branes; these strings do not feel the presence of the domain wall at all. Strings corresponding 
to other gauge fields, on the contrary, have at least one end on the fingerstall branes and 
are bound to them, so that they cannot propagate along the entire cylinder. On the other 
hand, we are not aware of the stringy interpretation of the localized fermion mode. 

We suggested a way to embed the (U(Ni) — U (A/2)) domain walls into M-theory. Namely, 
we suggested that they emerge as solutions of the matrix theory corresponding to the curved 
supergravity backgrounds which have the form of pp-waves with inhomogeneous three-form 
field strength. This effect can be thought of as a generalization of the Myers effect [27] . We 
considered two ways of introducing effects of curved backgrounds into the matrix model: one 
based on calculations at weak curvature, and the other relying on supersymmetry. These 
approaches lead to apparently different prescriptions in our case, however they both yield 
backgrounds admitting the domain wall solutions. The advantage of the latter approach is 
that it leads to the domain wall solution which saturates the BPS-type bound for the full 
matrix model action, and thus is stable. It is worth noting that our arguments rely on the 
approximation of weakly curved background; one may hope that they may apply beyond this 
approximation, especially taking into account that pp-wave backgrounds similar to those we 
discussed here were shown to be exact string backgrounds [22, 21]. We leave aside an issue 
of the supersymmetrization of the domain walls, though the BPS property suggests that it 
should be possible. 

To conclude, results obtained in this paper demonstrate that the rank of the gauge group 
can be a non-trivial dynamical parameter in the NC gauge theories. 
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